IV 


N90- 19944 

GRAVITATIONAL RADIATION FROM 
ROTATING GRAVITATIONAL COLLAPSE 

RICHARD F. STARK 
Institute for Theoretical Physics 
University of California, Santa Barbara, CA 93106, U.S.A 

abstract 





We find the efficiency of gravitational wave emission from 
axisymmetric rotating collapse to a black hole to be very low: AE I Me < 

7 x 1 0 -4 The main waveform shape is well defined and "^y 
independent of the details of the collapse. Such a signature will al low 
Dattern recognition techniques to be used when searching 
experimental data. These results (which can be scaled in mass) have 
been obtained using a fully general relativistic computer code that 
evolves rotating axisymmetric configurations and d y P 

their gravitational radiation emission. 

The results summarized here have come from a tested, fully general 

S-tHH ■£■ fS sS Sfes 

'a h „d p e i^%;^, C “4 JS. Son ““ nomaxisymmetric collapse 

is given in Stark 1988.) 

The axisvmmetric rotating gravitational collapse we have studied is that of an 

symmetric non-rotating polytrope. For the cases studied • h 

initial stellar radius of a 6 GAf c 2 and a fixed adiabatic index r = 2. We then reduce tne 
Xre -d . henna, energy to a fraction /, of their e q nil,br.nm values (w .h / = 

0 01 or 0 4) and simultaneously add a 'rigid body' azimuthal rotation to the star The 
rotation is measured by the dimensionless angular momentum p 

a - j I QM- L w here J is the total angular momentum of the star. The collapse and 

resulting gravitational wave emission are studied for a range of angular momenta 0 _ 
resulting ^avitaUor Qf state , a „ quantiti es scale in an elementary 

fashion wi.h M the mass of’ the star. Our initial conditions correspond to an tnmal 

radius of 8.8 x 10 5 M/M 0 cm and a central density of 1.9 x 10 (Ml M 0 ) & cm 
(corresponding to a somewhat unrealistic, compact, relativistic stellar core t ft* steMa 
masses). (For space-based experiments one can s.mply scale to h.gher M , 
approximating supermassive or population III collapse.) 

(For illustrations see the figures in Stark and Piran 1986). f.~ °^ h ^hJ 

f =0 4 ) the collapse takes place almost spherically symmetrically with the 

SZ ri:, * 

of a more rapidly rotating star with a - 0.75 (and f p ■ )> 


18 



become noticeable with the star soon becoming highly flattened into the equatorial 
plane. After a time of -30 M, the star bounces vertically away from the equator, but 
still continues to collapse radially inwards until freezing as a black hole forms, and 
the lapse collapses. For a > a cr i t (with a cr it = 0.85 :l 0.05 for f p - 0.4), the effects of 
rotation become dominant, no black hole forms, and the star ends up oscillating about 

a flattened rotating equilibrium configuration (see Stark and Piran 1986). (The 
general behavior for f p = 0.01 is similar to that described, except that the star 

becomes rather more flattened and (for a < a cr i t ) collapses slightly faster). 

Fig. 1 shows results obtained for the even and odd transverse traceless 

waveforms + (as defined in Stark and Piran 1987) from the collapse to black 

holes of stars with various a < a cr i t (and for f p = 0.01). The waveforms shown have 
been monitored at the outer edge of the grid (r - 50A/). The even and odd modes have 
the expected sin 2 0 and cos 9$in 2 6 angular dependences respectively, and these as well 
as the r -1 radial fall-off have been factored out. The waveform is well established by 
fairly small radii ( r - 25Af). The waveform consists of a broad peak at a retarded time 
-0 due to the initial flattening of the star from its initial spherical configuration. 

(This part contributes negligibly to the total energy of the gravitational wave 
emission.) The main emission appears at a retarded time of -25 M, corresponding to 
when the collapsing star is only ~3 M in size. The main emission lasts for -10A/ and is 
followed by a decaying oscillatory tail. The most noticeable feature of these 
waveforms is the insensitivity of the waveform shape to the value of a. Over the 

range 0 < a < a cr i t for black hole collapse, the shape remains closely unchanged, the 

waveform differing for each collapse only in its amplitude. (This is the case for f p = 

0.4 also.) 

The maximum h+ amplitude, | h+ \ max% from the collapse to a black hole, both for 
f p = 0.01 and 0.4, scales very closely as a 2 and finally levels off to a maximum value as 
a -> &crit • Our results are closely approximated by: 

(r IM ) | h+ I ma x = min{0.1a 2 , A max ) (0 < a < a cr i t ) 

with A max = 0.06, a crit = 1.2 ±0.2 for f p = 0.01; and A max = 0.025. a crit = 0.8 ± 0.05 for f p = 

0.4. (As would be expected for axisymmetric collapse | h x \ < 0.21/1+1.) 

The dependence of the total gravitational wave energy emitted, A E, on the 
angular momentum of the collapse (to a black hole) is shown in 

Fig. 2 for both f p = 0.01 and 0.4. The most noticeable feature is the very low efficiency 

cases the efficiency is < 7 x 10 -4 . The emitted energy A E scales as a 4 (as 

shown by the solid line fit in Fig. 2) with a coefficient of 1.4 x 10 -3 independent of f p , 
and A E levels off to a maximum value 6 max (which depends on f p , as shown by the 
dotted lines) as a -> a crit . Our results closely follow the form: 

A E I Me 2 = min{ 1 .4 x 10 _3 a 4 , 6 max } (0 < a < a cr it) 

with e max = 6 x 10 -4 iorf p = 0.01; and e max = 1 x 10 -4 for/j, = 0.4. In all cases (AE/Mc 2 ) + 

> 10(A ElMc\. 
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Corresponding to the insensitivity of the waveform shape, the shape of the 
energy spectra, from these collapses to a black hole is similarly insensitive. The 
spectra have the form of a simple peaked curve, although for values of a near a crit , a 
low frequency component due to the initial flattening of the star appears (see Stark 
and Piran 1986 for detailed figures). The energy spectra peak at frequencies of 

0.035 < v < 0.08 (^r) -1 Hz (i.e., 7 < v < 16 (M/Mq)’ 1 kHz) corresponds to wavelengths of 

between 12 and 28 ( 'M f) The peak spectral energy F(y max ) scales as a 4 , and levels off 
to a maximum value as a -> a crit . Our numerical results closely follow the relations: 

Fi^max) erg cm -2 Hz -1 = min{60a 4 , F max }(M/MQ) 2 (r/10kpc) -2 
with F max = 22 for f p = 0.01; and F max = 4 for f p = 0.4. 

Our waveforms are very similar in shape to those of several previous 
perturbation studies (see Stark and Piran 1986 for details). This broad agreement in 
waveform shapes suggests the importance of black hole quasi-normal mode 
excitation in the gravitational wave emission process. The collapsing star produces 
most gravitational wave emission when it is ~3 M in size, by which time the exterior 
metric is already quite similar to that of a black hole. The star thus acts as a source to 
excite the black hole quasi-normal modes. (The star being smaller than the peak of 
the black hole potential determining the normal modes.) The waveform shape 
reflects the properties of the black hole formed, whereas the amplitude to which the 
modes are excited depend on the details of the collapse. The observed insensitivity of 
the shape of the waveform to the rotation parameter a, corresponds to the weak 
dependence on a of the m = 0 axisymmetric modes of the Kerr black hole (Detweiler 
1980). This picture is confirmed by the very good fit we obtain to the waveform from 
a linear combination of the two lowest modes (see Stark and Piran 1986). 

We may summarize our results as follows. The most optimistic case ( i.e ., the 

one producing the largest gravitational radiation) is for a ~ a cr it an d a small initial 
pressure (f p = 0.01). For this case, we expect the maximum observable strain 
I A/// I max lo be: 


f \h+\ max = 4.0 x 10~ 19 (M/1.4M g ) (r/lOkpc) -1 
I Al/l\ max - |^/ g/c 3 F(Vmajc) = 3. 6 x iQ- 19 (Af/1.4M 0 ) (r/lOkpc)" 1 


where the two measures of the strain correspond to nonresonant (interferometer) 
and resonant (bar) experiments respectively. The maximum spectral energy, peak 
frequency, and full-width, half-maximum frequency are: 

F(v max ) = 45 (M/1.4M G ) 2 (r/10kpcr 2 erg cm -2 Hz -1 

v max = 10(A//1.4 Mq) - 1 kHz; Av J/2 = 4 (A//1.4M 0 ) _1 kHz 

For low a collapses ( a < 0.8), the expected maximum strain and maximum spectral 
energy are: 
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I A//Z \ max = 6.7 x 10" 19 a 2 (A//1.4M G )(r/10kpcr 1 


(both cases) 


F(v m ax)= 120 a 4 (M/l .4M G ) 2 (r/10kpc) 2 erg cm' 2 Hz' 1 
with v max and A v 1/2 as above. 

An important feature of the main black hole emission is that the general 
waveform shape is well defined and nearly independent of the details of the collapse. 
This provides a signature for the black hole collapse and will allow pattern 
recognition techniques to be applied to broadband experimental gravitational wave 
data, thus allowing greater sensitivities to be reached for these events. Further 

experimental implications for ground-based bar detectors of stellar collapse can be 
found in Stark and Piran 1986. 
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FIG. 1 — /*+ (left) and hx (right) waveforms vs. retarded time for various a (/ = 0.01) black hole 
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FIG. 2 — Total Gravitational wave energy (ADM) emitted vs. a for f p = 0.01 (squares) and / 
(triangles) black hole collapses. Solid line fits the a 4 scaling; dashed lines mark e 






DISCUSSION 


ARMSTRONG: How would you expect the shape of gravitational temporal waveforms to 

change for non-axisymmetric collapse? Gravity wave efficiencies? 

STARK: The waveforms from rotating non-axisymmetric black hole collapse can be 

expected to consist of a quasi-periodic phase (due to the rotation of the bar) followed 
by emission from the actual black hole formation and a decaying tail. Th.e ma ' n 
black hole emission can be expected to show similar normal mode characteristics as 
we have seen in the axisymmetric collapse. Thus we may expect this part of the 
emission to have a shape similar to that for the axisymmetric collapse, except possible 
at high rotations where the non-axisymmetric normal mode frequencies become 

appreciably split. As far as the efficiency is concerned, perturbation studies indicate 

that it may be at least an order of magnitude larger than for axisymmetric collapse 
(for an appreciable level of non-axisymmetry). We will not know for sure until 

numerical simulations are made. (See e.g., Stark, 1988 for work in progress). 

HELLINGS: Is it true that changing "a" changes the amplitude without changing 

the period while changing ”m" changes both in a proportional way? 

STARK: The waveform amplitude from collapse to a black hole scales as a 4 for low a 

(where a = J/GM 2 /c; J the total angular momentum; M the mass of the star) and then 

levels off before a reaches a C rit (the rotation above which black hole formation no 
longer occurs). The amplitude also scales proportionally with M. The waveform 

period (or more precisely, the scale of retarded time for the waveform) of the main 

emission from black hole formation is insensitive to a (due to the insensitivity of the 
axisymmetric Kerr black hole normal modes frequencies to a). The period scales 
proportionally with M. 

SHAPIRO: Given the complexity of the system you're modelling, and the (necessarily) 

complicated code needed to solve the coupled equations used to represent the system's 
evolution, what checks have been made to try to insure that the results are correct? 

STARK: It is clearly very important that a code of this kind be tested as fully as 

possible. We have performed an extensive series of tests which can be broadly 

classified into general stability tests, conservation tests, hydrodynamic tests and 
comparison tests against known approximate perturbation solutions. (The derivation 

of the equations we solve has also been checked using symbolic computer algebra). 

A detailed description of these tests can be found in Stark & Piran, 1987. The tests 

performed include: (i) Evolution of the vacuum Schwarzschild exterior (r > 2M) for 

many gravitational times; (ii) Propagation of gravitational waves (for both 

polarizations) along inward and outward characteristics with negligible reflection at 
the outer boundary; evolution to flat spacetime for low amplitude waves, and black 
hole formation for high initial wave amplitudes; (iii) Stable evolution of initial data 

consisting of flat spacetime + a small amplitude 'random data’ superimposed; (iv) 

Evolution of stable extreme relativistic polytropes for many free fall timescales; 

including checking the excitation of the lowest radial mode oscillations for adiabatic 

indices G > G C rit< as well as the instability to collapse for G < G cr i t (with G cr it = 4/3 + 
general relativistic corrections); (v) Conservation of the ADM mass and the total 

angular momentum; (vi) Conservation of the specific angular momentum spectrum, 
(vii) Propagation of generalized linearized Teukolsky waves for both polarization 
modes; (viii) Comparison of the gravitational wave emission from the infall of a 
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spheroidal dust shell onto a Schwarzschild black hole with known perturbation 
results, and (ix) Comparison of the gravitational wave emission and hydrodynamics 
lor uniformly rotating homogeneous spheroidal collapse with known Newtonian + 

fr^^ 0 C r ^? rma ! Sm r , esi ! Its ' 11 ,s also worth remarking that the waveform obtained 
the full scale calculation of the axisymmetric collapse to a black hole shows a 
remarkable agreement with the waveform from an approximate perturbation study 
(Nakamura Oohara & Kojima, Prog. Theor. Phys. Suppl., No. 90, 1987) of a rotating 
dust ring falling in from infinity onto a Schwarzschild black hole with a specific 
angular momentum. comparable to that of the collapsing star. The collapsing dust 
ring mimics the flattening in the equatorial plane of the collapsing star. One finds 

very near agreement in the waveform shapes, and, on scaling up the perturbation 
result (beyond its actual range of validity) the amplitudes as well. (Both efficiencies 
then also agree at a few x 10' 4 ). 

BENDER: How likely does it seem that some instability such as a bar instability would 

occur during collapse? What mechanisms might determine the initial axial 

asymmetry ? 

STARK: A bar mode instability is very likely for sufficiently high rotation The 

amount of non-axisymmetry depends on the exact details of the microphysics and 
initial data. At present these are not known sufficiently well. Work in progress on 

non-axisymmetric collapse to a black hole (Stark, 1988) will investigate how the 
gravitational emission varies for a range of non-axisymmetry, but will not attempt to 
answer which level of non-axisymmetry is astrophysically appropriate. 

SCHUTZ. I will just comment on the question of the generation of non-axisymmetry 
m collapse Because it relies on an instability that grows exponentially, only a very 

small initial perturbation is needed. The instability is the ’bar mode' instability that 

has been discussed as well in the context of formation of bars and spirals in galaxies 
Basically, if angular momentum is conserved in the collapse it becomes 
energetically favorable for the system to deform into a tumbling^. The hard 
question for numerical calculations to answer, and for which careful attention to the 
microphysics may be necessary, is how much angular momentum transport there 
may be as the object collapses. K 
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